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We calculate the two-point correlation function using the interpolating current with J^ = 2 . 
After performing the Borel sum rule analysis, the extracted masses of the 2^^ tensor charmonium 
and bottomonium are 3.89 ± 0.08 GeV and 10.06 ± 0.09 GeV respectively. For coinparison, we 
also perform the moment sum rule analysis for the charmonium and bottomonium systems. As a 
byproduct, we extend the same analysis to study the qq,qs,ss,qc,sc,qb,sb and cb systems. Their 
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I. INTRODUCTION 



(N 

JL , Charmonium spectroscopy provides a crucial test of the quantum chromodynaniics(QCD) in both the perturbative 

Qj ' and nonperturbative regimes. The charmonium spectrum can be calculated using the potential models [l|. In the 

r^ [ picture of the conventional quark model, the charmonium states are characterized by the J^ quantum numbers: 

P — (— )^+^,C = {~)^~^ : where L is the orbital angular momentum and S the total spin. Their quantum numbers 

are J^^ = 0±+, 1^", 1++, 2++, 2-± and so on. 

^ , There has been a revival of the charmonium spectroscopy in the past eight years. Since the operation of the large 

^s^ ' facilities such as Tevatron and two B-factories, many unexpected charmonium or charmonium-like states above the 

T^ ■ open charm threshold have been discovered, such as X(3872), r(4260), r(4360), r(4660), Z+(4430) etc. [i-[|. Some 

0^ ' of them tend to decay into a charmonium state plus light hadrons. The underlying structures of these new states are 

'^ ] not known precisely. They are sometimes speculated to be the candidates of the exotic states such as the molecular 

f — ■ states, the tetraquark states, the hybrid charmonium, the baryoniuni states and so on. 

^^ ■ After the discoveries of rjc{2S) and hc{lP) 0, Qi all the charmonium states below the open charm threshold have 

been observed experimentally @. They are J/^{l^Si),i:{2S){2^Si),_Xco,i,2{l^Po,i,2),Vc{'i-^So),Vc{2S){2^So) and 
hc{i^Pi)- All these charmonium resonances are narrow. Below the DD threshold, there is only one tensor meson 

Xc2(lP). 

In the study of the meson spectroscopy, the local meson interpolating fields iJj{x)T'4>(x) are usually introduced. 
. , However, these operators are useful only in the study of the low-lying states with J^'^ = O^"*", 1^~, 1+"'" etc. In 
jrt ' order to explore the higher spin states, the non-local fermion operators with covariant derivatives acting on the quark 
fields should be used. For example, the authors studied the charmonium spectrum including the higher spin states 
by including the non-local operators in the framework of the lattice QCD simulations in Ref. [lO|. The masses and 
decay constants of the ground states heavy XQ2 tensor mesons have been calculated by using the QCD sum rules 
approach in Ref. [ll|. The newly observed resonance X(1600) was studied as a J^'^I'^ ~ 2++2"'" four- quar k state in 
the framework of QCD finite energy sum rules in Ref. [12]. The tensor currents with J^'~^ — 2++ [ij, [iJl and 2~~ 
141 were firstly employed to study the light quarks systems. 

In this work, wc use the tensor current with J^^ = 2~~ to calculate the two-point correlation function in the 
most general situation. Then we perform the QCD sum rule analysis and extract the masses of the 2 tensor 
states. Especially for the charmonium and bottomonium states, we also perform the moment sum rule analysis for 
comparison. None of these possible resonances has been observed up to now except for the strange mesons iir2(1770) 
and i<r2(1820), which have I{J^) = 1(2") with no definite C parity. 
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The paper is organized as follows. In Sec. |TT1 we discuss the quantum numbers of the interpolating current and 
calculate the two-point correlation function in the general situation. In order to cross-check the quantum number of 
the interpolating current, we discuss the reduction of the current in the non-relativistic limit in Appendix [A] We 
perform the Borel sum rule analysis for various systems in Sec. IIIII For comparison, we present the moment sum rule 
analysis for the charmonium and bottomonium systems in Appendix [C] The last section is a short summary. 

II. THE TWO-POINT CORRELATION FUNCTION 

In the framework of QCD sum rule [l5|-[l7|, we consider the two-point correlation function: 

n^.,p.(g') = I J d''xe^'^^^~y\0\T[J^,{x)jl{y)]\0) 
1 2 

= -^{vi^pV^a + i^'i'nyp - 2^M>'^p'^)n('7^) + •■•, (1) 

with rj^i, — qi_iqu/(f' — g^u- Jiiv is the tensor interpolating current with J^*^ = 2 : 

Jpu = Ql{x){"^^i-i5Du + lul5uli - ^Vt^i^l5lJ)Q2ix). (2) 

The rjfj^i, term is introduced to ensure the trace condition: 

5^" V - 0' (3) 



The covariant derivative _D^ is defined as: 



^^ = Dt-K, (4) 



D, 



d', + zg-AlX = K-^9^Al (5) 



The interpolating current in Eq. ^ was first constructed in Ref. [lj| to study the tensor meson /2(1670) in the 
framework of QCD sum rule. It was also introduced in Ref. [l^ as an operator with J^^ = 2 ^ . However, we will 
show in Appendix!^ that it carries the odd C parity through the charge conjugation transformation. Wc also discuss 
this tensor operator in the framework of the quark model. In fact, its J^ component reduces to the D-wave tensor 
operator while its Joi component reduces to the P-wave axial- vector operator in the non-relativistic limit. 

In order to study the tensor resonance, one should pick out the intrinsic spin-2 tensor structure from the two-point 
correlation function induced by the tensor current. In Eq. ([T]), n(q^) contains contributions from the pure tensor states 
only and "..." represents the other possible structures from other states such as the spin-1 states. At the hadron level, 
the correlation function Tl{q^) satisfies the dispersion relation: 

n(.^) = r ^^, (6) 

The hadronic spectral density p(s) is usually assumed to take the pole plus continuum contribution: 

p{s) = Y,S{s-ml){QMn){n\rj+\0) 

n 

~ /3J-TO3f(5(s — 771y) + continuum, (7) 

where mx denotes the mass of the resonance X and fx stands for the coupling constant of the tensor meson to the 
current J^j^: 

(0|V|X)=/xm|e^, . (8) 

Here e^^ is the polarization tensor. 

The correlation function can also be computed at the quark-gluon level using the operator product expansion(OPE) 
method. For the heavy quark systems, it's convenient to evaluate the Wilson coefficient in the momentum space. 
In our calculation we only consider the first order perturbative and various condensates contributions, i.e., the bare 



quark loop in Fig.[T] {g'^G'^) terms in Fig.[2]and (g^fG^) terms in Fig.|3l The massive quark propagator in an external 
field Afj^ (x) in the fixed-point gauge is listed in Appendix |B] The quark lines attached with gluon legs contain terms 
proportional to y. We keep these terms throughout the evaluation and let y go to zero after finishing the derivatives. 
In Fig. [5] and Fig. [31 the diagrams with a gluon leg attached at the right vertex are linear in y and vanish after putting 
y = 0. They do not contribute to the correlation function. 




FIG. 1: Feynman diagrams of the perturbative term for the vacuum polarization, 
line and interpolating current respectively. 



The solid and dashed lines denote the quark 




FIG. 2: Feynman diagrams of the {glC^) contribution to the vacuum polarization. The solid, curly and dashed lines denote 
the quark, ghion lines and currents respectively. 




FIG. 3: Feynman diagrams of the (pJ/G^) contribution to the vacuum polarization. 












FIG. 4: Feynman diagrams of the (gsqq)^ contribution to the vacuum polarization. 



In order to suppress the higher states contributions, it is significant to perform the Borel transformation to the 
correlation function, which also helps improve the convergence of the OPE series. With the assumption of the 
quark-hadron duality, we derive the tensor meson sum rule: 



n(M^) = /im^e 



6^ -ml /Ml 



dser'/^^Bpi^s), 



' (??ii+rn2)^ 

where sg is the threshold parameter. Then we can extract the meson mass mx'- 



(9) 



mx = -- 



d 

a{i/M%[ 



■n(A4) i;i^....^.dse~^/^'lsp{s) 



+m2)^ 



n(A/|) 






After performing the Borel transformation, the correlation function reads: 
_6 

"" J(mi+m2) 



(10) 



np'=''*(M|) 



if 



nr '{M 



n 



{G^> 



{Ml) = 



247r2 
247r2 



(mi+m2) 
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dx 





dse-^'^'', 
(4x3 



dx 



< [rn^ — s(l — x)] X + 01,2(1 — x) > [77111712 — {m\ — mj + sx)x\ , 



2x' 



dx{Sx - 4x - 5), 
9x + 2)t7iI 2(9x2 - 9a; + 4)mim2 (4^^ - 4x - 5)m 



(1-x) 



(1 — x)x 



2x(x + 2)r7i^ {2x + 3)777^771-2 2{x:^ — x + l)m\m 
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M|(l 
2to| 
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M2x2 



M|x2 



n<«'>(M|) 



2887r2 



dx 



(20x3 + 3x2 _ gg^ , 



35) 3(2x2 ^ 3^ 



2)TOf 2(5x2 + 4x-2)mfm2 



(l-x)2 M|(l 

(10x3 - 32x2 + 39x - 14)mfm 



+ 



6(x2 



Ml (1 - x)3x2 

4x + 2)mi'ml (32x'^ 



4 M|(l - x)3x 

(32x4 _ 7;^^3 _^ ^4^.2 + 6x - 7)to? 

A/|(l-x)3x 



46x3 + 124x2 _ gg^ ^ 33)mim2 



^1 

3?7li 
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(32x3 - 87x2 _^ 72x - 20)m: 
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B 

-ml/M 



Ml{l-x)'^x 



2m,i{q2gsO- ■ Gq2)e 



By 

-val/Ml 



'/Ml 



(11) 



where Xmin = i <! 1 

m^ = ^ ^^ ^ 



i|l + IIiI_I2i_ [(l + Z2l^)2_4„,2/,]i/2|^^^_^ ^ i|l^™i^ 



and 



,J^_^^ . Only the two-gluon condensate contribution H^'^ H^^s) and tri-gluon condensate contribution 

n'"^ ' [Ml) are involved for the heavy quark systems (Qi, Q2 = c, 6). For the gg, qs, ss, qc, sc, qb and sb systems, the 
four quark condensate g1{qiqi){q2(l2) (only for the light quarks systems, (71,92 = u,d,s) and the quark-gluon mixed 
condensate {qgs<y ■ Gq) arc also needed. Among these systems, the charge neutral ones have the quantum numbers 
J^'-^ = 2^^ and the charged ones have J^ = 2^ with no definite C parity. The quark condensate {qq) does not 
contribute to the intrinsic tensor structure. The four quark condensate g1{qiqi){q2(l2) as shown in Fig. |4l always 
plays an important role in the conventional qq meson sum rules |15l. Il6|. In the present case, only the first diagram 
in Fig. |4] contributes to the correlation function. All the other diagrams vanish due to the special Lorentz structure 
of the current. 



III. NUMERICAL ANALYSIS 



In the QCD sum rule analysis we use the following values of the quark masses and various condensates |9l. ITsI, ITg. [Tsj : 

TO„(2 GeV) = (2.9 ± 0.6) MeV , 
md(2 GeV) = (5.2 ± 0.9) MeV , 
TO, (2 GeV) = (4.0 ± 0.7) MeV , 
TOs(2 GcV) = (lOlla?) MeV, 
{qq) =: -(0.23 ±0.03)^ GeV^ 
{qgsOGq) = ~Ml{qq) , 



M^ = (0.8 ± 0.2) GeV^ , 
{rs.s)/{qq)^0.8±0.2, 

{glGG) = 0.88 GeV^ , 
{glfGGG) = 0.087 GcV^ 



(12) 



where the up, down and strange quark masses are the current quark masses in a mass-independent subtraction scheme 
such as MS at a scale /z = 2 GeV. The running charm quark mass has been determined by the moment sum rule 
in Refs. Il9l-j23l |. Recently, the value has been updated by using the four loop results for the vacuum polarization 
function |24l |. 

Since we have not calculated the next leading order radiative correction due to the complicated interpolating current 
in the present work, it is desirable to extract the charm quark mass within the same QCD sum rule formalism using 
the experimental J /ip mass as input. The J /^ sum rule derived from the interpolating current j^ = C7^c was known 
very well [la ] . With the same criterion of the present tensor current and keeping only the leading order perturbative 
term, gluon condensate and tri-gluon condensate contributions, we show the Borel sum rule results in Fig. [5] using 
the experimental data [9|. The extracted charm and bottom quark mass are nic = 1.35 GeV and mi, = 4.60 GeV as 
shown in Fig. [5] 
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FIG. 5: The extraction of the running charm and bottom quark masses in Borel sum rules. 



After performing the Borel transformation, the correlation function in Eq. (fTTj) is the function of the threshold 
parameter sq and Borel mass Mb- In the Borel sum rules analysis, there should exist suitable working regions of 
these two parameters in order to obtain a stable mass sum rules. For this purpose we choose the value of sq around 
which the extracted mass mx is stable with the variation of Mg. In Eq. pTj) . the exponential weight functions 
suppress the higher states contributions for the small value of Mg. However, the convergence of the OPE series 
becomes bad if Mg is too small. These two opposite requirements restrict the domain of the Borel mass. We define 
the pole contribution (PC) as: 






The constraint of a relatively large pole contribution leads to the upper bound M, 
the requirement of the OPE convergence yields the lower bound M^^^. 



2 
max 



(13) 



of the Borel parameter while 



Using the parameter values in Eq. p^ . we study the tensor charmonium system by considering only nf"^''*(Af|), 
n^'^ )(Af^) and 11^'-' \M]j) in Eq. (fTT]) with toi = m2 = rric. The two-gluon condensate contribution is the dominant 
correction to the correlation function in this situation. The lower bound of the Borel parameter is obtained as 
■^min — 2-5 GeV^ by requiring that the two-gluon condensate correction be less than one fifth of the perturbative 
term and the tri-gluon condensate correction less than one fifth of the two-gluon condensate correction. The upper 
limit of M^ is the function of Sq as shown in Eq. (|13p . We choose ^/sq = 4.3 GeV around which the variation of the 
extracted mass mx with Mg is the minimum, as shown in Fig. |B^. Then we require that the pole contribution be 
larger than 40% to get the upper bound of the Borel mass M^^^^ = 4.3 GeV^. 

By performing the numerical analysis in the domain 2.5 < Mg < 4.3 GeV^, we obtain a stable mass sum rule of the 
2 charmonium system. The dependence of the extracted mass mx with the Borel parameter is very weak in this 
domain of Af^, as shown in Fig.[B}D. The extracted mass for the possible charmonium tensor state is about 3.89 ±0.08 
GeV. As expected naively in the quark model, this value is slightly higher than the mass of the lowest D-wavc 1 
charmonium state "0(3770). 

We extend the same analysis to the cb and bb heavy quark systems and collect the numerical results in Table HI The 
masses of the Be and bottomonium tensor states are extracted to be 7.01 ± 0.09 and 10.06 ± 0.09 GeV, respectively. 
Only the errors from the uncertainty of the threshold values and variation of the Borel parameter are taken into 
account. Other possible error sources include the truncation of the OPE series and the uncertainty of the quark 
masses, condensate values and so on. The mass of the 2~ cb state was predicted to be around 7.0 — 7.1 GeV in 
Refs. |25l - [28| . which is consistent with our calculation. 
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FIG. 6: The variation of mx with so(a) and Mg{h) for the 2 charmonium state in the framework of the Borel sum rules. 
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FIG. 10: The variation of rax with sq and M% for the 2 qc state in the framework of the Borel sum rules. 
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FIG. 11: The variation of rax with so and M% for the 2 sc state in the framework of the Borel sum rules. 
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FIG. 12: The variation of mx with so and Mg for the 2 qb state in the framework of the Borel sum rules. 
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FIG. 13: The variation of mx with so and Mg for the 2 sb state in the framework of the Borel sum rules. 
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FIG. 14: The variation of rax with so and M% for the 2 ch state in the framework of the Borel sum rules. 
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FIG. 15: The variation of mx with so and M% for the 2 bottomonium state in the framework of the Borel sum rules. 



We also extend the analysis to the qq, qs, ss, qc, sc, qb and sb systems, where q denotes u or d quark. As mentioned 
in Sec. [Ill the four quark condensate g1{qiqi) {q2q2) and the quark-gluon mixed condensate (ggsC ■ Gg) should be 
considered now. The corresponding parameters such as the quark masses and the condensates should be used for 
different systems in Eq. (jlip . Especially for the qc {sc) system, only the light (strange) quark-gluon mixed condensate 
needs to be considered. After performing the numerical analysis, the variations of the mass with the threshold value sq 
and Borel parameter Af^ are shown in Figs. IMTSl We show the Borel window, the threshold value and the extracted 
mass for all these systems in Table IH As shown in Table HI the light tensor meson mass is about 1.71 GcV, which is 
consistent with the value 1.63 GcV in Ref. \\M. 
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^4.3 


3.89 ±0.08 


qb 


6.22 35^ 


^5.3 


5.72 ±0.11 


sb 


6.92 4.0- 


-8.4 


6.37 ±0.10 


cb 


7.42 4.0 - 


-7.0 


7.01 ±0.09 


bb 


10.42 5.0- 


-8.5 


10.06 ±0.09 



TABLE I: The threshold value, Borel window and extracted mass for the various 2 tensor states in the framework of the Borel 
sum rules. 

IV. SUMMARY 

We have shown that the interpolating tensor current in Eq. ^ is the correct operator for the J^^ = 2 meson, the 
Jij and Joi components of which reduce to the D-wavc and P-wave states respectively in the non-relativistic limit in 
terms of the quark model picture as shown in the appendix. Then we make the operator product expansion (OPE) and 
calculate the two-point correlation function. For the heavy quark Q1Q2 (Qi, Q2 = c, b) systems, the power corrections 
include the gluon condensate {giGG) and the tri-gluon condensate (g^fGGG) only. For the light quarks systems and 
the light-heavy quarks systems, the four quark condensate g1{qiqi){q2q2) (only for the light quarks systems) and the 
quark-gluon mixed condensate {qgs<J ■ Gq) also contribute to the sum rules. While the gluon condensate (g^GG) is 
the dominant nonpcrturbative correction for the above systems, these terms also play an important role in the tensor 
meson sum rules. 

Within the framework of the Borel sum rules, we have studied the qq,qs,qc,sc,cc,qb,sb,cb and bb systems. All 
these systems display stable QCD sum rules in the working regions of the Borel parameter Mg. Up to now, none 
of these 2~ tensor mesons have been observed except the strange meson iir2(1770) and 7^2(1820) Ql- As shown 
in Table. HI the extracted mass of the qs tensor state is about 1.78 GeV, which is consistent with the mass of the 
strange mesons 1^2(1770) and ii:2(1820) d]. The lowest D-wave charmonium state is ^^(3770) with J^^ = 1 — . The 
extracted mass of the 2 D-wave charmonium state is around 3.89 GeV, which is slightly higher than ^'(3770) as 
expected in the quark model. For the charmonium and bottomonium systems, we also perform the moment sum rules 
analysis for comparison. Hopefully the present investigation will be helpful to the future experimental search of these 
tensor states. 
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Appendix A: THE QUANTUM NUMBERS OF THE INTERPOLATING CURRENT 

In order to study the quantum numbers of the tensor current in Eq. ([2]), we perform the parity transformation and 
charge conjugation transformation to Jf^^: 

C^C-i = -V (Al) 

where (—1)^ = 1 for /x = and (—1)^ = —1 for ji = 1,2,3. With these relations and the trace condition in Eq. ([3|), 
the tensor current can couple to the 1^ (Joi components, i=l,2,3) and 2 (J^ components) states at the same 
time. 

It's interesting to reduce this operator in the non-relativistic limit and center of mass frame in terms of the quark 
model language. 

_ ^ Y 

Jij{x) = Qi(a;)7i75 dj Q2{x) + (i o j) 



i^i-^iiH^) (1. ^) f? i) f^^i i^kM^^j) 



Xr 



u{p' , s)'-fij5ikjv{p, r) + (i O j) 

'Ep,+m' \~'J^ Oy) 1^1 

■^(filicTtakPk + crkcriPk)Xr + {i ^ j) 

17^V'l[('^ifc + i<^ikim)Pk + {hi + i<^kiicri)Pk]Xr + (« ^ j) 

— ipl[k.jqi + ieiki(Jikkkj]xr + (* ^ j) 



--T—(pleikicrikkkjXr + {i^ j) (A2) 

Zm 



11 



Joi{x) 






^ t a-p' , (Q 1 






Xr 



(iki 



(A3) 



Joo -- ^fcoV'lXr = 



(A4) 



in which k ~ p — p' and q = p + jf = in the center of mass system. We have used the non-rclativistic hmit: 
Epi = Ep = m,ko = 0. It is obvious that Jij reduces to the D-wave and Joi reduces to P-wave in the non-relativistic 
hmit. Therefor the quantum numbers of the current should be 2 for the Jij component and 1"' for the Jqi 
component. 

Appendix B: THE MOMENTUM SPACE PROPAGATOR 

To calculate the higher dimensional giuonic operators, we consider the giuon field as an external one with the 
fixed-point gauge condition |29l - [3l| : 



ix-xorAl{x)^0, 



(Bl) 



where xq is an arbitrary point in space which can be chosen to be the origin. Then the four potential A" can be 
expressed in terms of the field strength tensor Gn,y{G^i, = ^G" ): 



M^) 



tdtGi,fj_{tx)x'^ 



Ix-'G^^iO) + ^x^x^D^G^^iO) + ^x^x^x'^D^DpG.^iO) 



» 
o 
o 
o 
o 
o 
o 



(a) 



(W 





FIG. 16: The graphic representations of the massive quark propagators with gluon legs attached to the quark line. 

Denote the massive quark propagator between the position x and y in the coordinate space as iSo{x, y). The massive 
quark propagator in the momentum space can be obtained as: 



iS{p) = iSoip) + iSg{p) + iSgg{p) + iSgggip), 

where iSo{p) is the free quark propagator as shown in Fig. I16r a): 

i 



iSo{p) 



p — m 



iSg(j>) is the quark propagator with one gluon leg attached as shown in Fig |16r b): 



. zA" ^„ a^^-(p + m) + ip + m)a^^'^ , z A" 



4 2 



(p2 — m^) 



2 2 



2y'^p''{p + m) y^^Y 



(p^ — m?) 



2^2 



p^ — m? 



iSgg{p) is the quark propagator with two gluon legs attached as shown in Fig 116( c): 



^Sggip) = ' — ^giGl^pG 



'lip^va (2 



4 2 2 



2r<a f.b P + ^ 

a r~ih 

l^P^y-y (p2 _ ^2)4 



[jp- — m?Y 
p + m 



t-ctippu _|_ Mpvpa j^ jpy<Tp\ 



y^r" + F""") + y'r""" - ly'y^f^^p' - m') 



(B2) 
(B3) 

(B4) 



(B5) 
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iSggg{p) IS tlic quRrk propagator with three gluon legs attached as shown in FigfTBld): 



iSgggip) 



i X" X'' A= 



p + m 






+F 



fiuPp-ya I j^fii^aplS-y i ^puap^fi i j^pvpa^-y , rpiyppa-y , j^pi^p^-ya. i j^fiupajP _, j^pup-yap 



F 



+ F 



i \'^ \b \<:- 

I A A A 3^a r<i> r"^ 



+ / 

p + m 



r 



+ F 



+ r 



+ r 



lyOi(ftJ-t^PPl _|_ jiJ-vpPi _|_ jp.upip\ _^^yPr-e 



13/ rp.avp-y 



rp.uap^ 



2 2 2"""^""'^'" ^^(^2-^2)6 

,rp.upja , rfj,upaj\ I j^yl/jP-avfip , rfiaupfl , rfj,uaPp , rpujlap . rfj.ul3pa , ^puapP , rfj,upaP , rpupPa\ 

-ly'^y^yT'^ip^^mY] (B6) 

where /A'''-"/^ = ^^(^2 _ TO2)^,.(-p2 _ ^2y^^^a(^p2 _ ^2)^/3(p2 _ ^2y 



Appendix C: Moment sum rule analysis 

For comparison, wc may also use the method of the moment sum rule |15l . |l6[ to study the 2 charmonium and 
bottomonium systems. To suppress the contribution of the higher states and pick out the lowest lying resonance, we 
define the moment by taking derivatives of the polarization function Il{q^) in Euclidean region q^ = —Q^ < 0: 



-i2n _ 1 / C? \ TJfr^2^ 



"•■<«' -;3 1 -5^J"«''i«--o-£i7Ti^ 



ds. 



With the spectral function in Eq. [71 we obtain: 

MniQl) -- 



fx^% 



{m\ + Ql) 



2\n+ 



j[l + 5^[Ql)l 



(CI) 



(C2) 



where i5„(Q§) denotes the higher states contributions to the moment divided by the lowest lying resonance contribution. 
To eliminate fx in Eq. (jC2p . one can consider the ratio of the moments: 



r{n,Qo) = 



MniQl 






Mn+l{Ql) '"^^ ' ""'l+<5„+i(Q2 



(C3) 



From the ratio r(n, Qq) one can immediately extract the mass of the lowest lying resonance nix for Sn{Qo) — ^n+i{Qo) 
at sufficiently high n. 

However, taking higher derivative n means moving further away from the asymptotically free region. This can be 
compensated by choosing a larger Qg. In Eq. (|C3p . it will be difficult to extract the mass of the lowest lying resonance 
for large value of Qq because (Sn(Qo) converges less fast in this situation. In fact, one can arrive a region in the (n, Qq) 
plane where the lowest lying resonance dominates the integral Eq. (|C1[) and the nonperturbative contribution is not 
too large. The moment Mn{Q^) can be drawn from the Borel transformed correlation function shown in Sec. HIl after 
taking into account the gluon condensate and tri-gluon condensate. 
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FIG. 17: The extraction of the running charm and bottom quark masses in moment sum rules. 
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Using the parameters in Eq. ([12)) and the experimental J/tp and T(IS') masses as input, we extraet the eharm and 
bottom quark masses within the moment sum rule framework as shown in Fig. 1171 m^. = 1.29 GeV and tti^ = 4.55 
GcV. Then we perform numerical analysis to extract the mass mx as the function of Qq and n. In Fig. 1181 we 
show the masses of the 2 charmonium and bottomonium states as the function of n. There is a mass minimum 
in these curves under the variation with n. By choosing Q^ = 50?7i^,60to^, 70to^ for the charmonium system and 



Ql = 



4m2,i 



'b' 



12to^ for the bottomonium system, the extracted mass converges to 4.03 GeV and 10.16 GeV 



respectively. These values are consistent with the results of the Borel sum rules within the errors. 
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FIG. 18: The variation of mx with n for the 2 
sum rules. 
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